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Abstract: We investigate the decomposability of nonnegative com¬ 
pact r-potent operators on a separable Hilbert space C 2 (X). We 
provide a constructive algorithm to prove that basis functions of 
range spaces of nonnegative r-potent operators can be chosen to be 
all nonnegative and mutually orthogonal. We use this orthogonal¬ 
ity to establish that nonnegative compact r-potent operators with 
range spaces of dimension strictly greater than r — 1 are decompos¬ 
able. 
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1. Introduction 

An operator A is said to be idempotent if it satisfies A 2 = A. 
An operator A on the Hilbert space C 2 (X) is said to be decompos¬ 
able DE] if there exists a nontrivial standard subspace of C 2 (X) 
that is invariant under A. Marwaha [3] showed that nonnegative 
idempotent operators of rank greater than one are decomposable in 
finite dimensions. Marwaha jUJ further established that nonnega¬ 
tive idempotent operators with range spaces of dimension greater 
than one are decomposable in infinite dimensions. The results in [3j 
(finite-dimension case) for idempotent operators were generalised 
by Thukral and Marwaha to r-potent operators in [5] (recall that 
an operator A is said to be r-potent [6] if A r = A, where r is a posi¬ 
tive integer). Specifically, Thukral and Marwaha [5] showed that all 
nonnegative r-potent operators of rank greater than r—1 are decom¬ 
posable in finite dimensions. A similar generalisation from idempo¬ 
tent operators to r-potent operators in the infinite-dimension case 
[4j however does not exist in the literature. In this paper, our 
main contribution is to provide this generalisation of results in [4] 
for idempotent operators to the r-potent case. It should be noted 
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that proof of decomposability for r-potent operators in the finite- 
dimension case [5] depends critically on the Perron-Frobenius The¬ 
orem [7] for matrices. Since there is no direct analogue of Perron- 
Frobenius theorem in the infinite-dimension case, generalisation of 
results in [3] to the r-potent case is completely nontrivial. In fact, a 
substantial portion of this paper is devoted to developing new tools 
and techniques that are eventually deployed to prove the decom¬ 
posability of r-potent operators in the infinite-dimension case. The 
rest of this paper is organised as follows. Section 0] contains the 
notations, definitions and preliminaries used in this paper. In Sec¬ 
tion El we provide a novel constructive algorithm to show that the 
basis functions of the range space of a nonnegative r-potent oper¬ 
ator can be chosen to be all nonnegative and mutually orthogonal. 
Finally, Section 0] states the result on decomposability of r-potent 
operators and the corresponding proofs. 

2. Preliminaries 

In this paper, X denotes a separable, locally compact Hausdorff 
space and /j a Borel measure on X. C 2 (X) is the Hilbert space 
of (equivalence classes) complex-valued measurable and square in- 
tegrable functions on X. For sake of simplicity, we assume that 
fJ,(X) < oo. 

Definition 1. [HI p. 57] A function / e C 2 (X) is said to be 
nonnegative, written as / ^ 0, if p{x e X : f(x) < 0} = 0. 

Definition 2. [12] A standard subspace of C 2 (X) is a norm-closed 
linear manifold in C 2 {X) of the form 

C 2 (U) = {/ 6 C 2 (X) : / = 0 a.e. on U c } 

for some Borel subset U of X. This space is nontrivial if p{U) ■ 
p{U c ) > 0. 

Definition 3. |2] An operator A on C 2 (X) is said to be de¬ 

composable if there exists a nontrivial standard subspace of C 2 (X) 
invariant under A. 

In this paper, we shall use the following equivalent definition of 
decomposability (provided as a proposition in [3J p. 39]): 

Definition 4. A nonnegative operator A on C 2 (X) is decom¬ 
posable if and only if there exists a Borel subset U of X with 
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g{U) ■ g{U c ) > 0 such that (A xu, Xu c ) = 0, where Xu = 1 on 
Supp U and Xu = 0 on Supp U c (Note that U c denotes the com¬ 
plement of U). 

Definition 5. Suppose X\ and X 2 are Borel subsets of X. An 
operator A from C 2 (X i) to jC 2 (X 2 ) is called nonnegative if A/ ^ 0 
whenever 0 ^ ^ 0 in £ 2 (dr). Similarly, A is called positive if 

A/ > 0 whenever 0 * / ^ 0 in C 2 (X 1 ). 

Definition 6. For any function /, we define the support of / 
as Supp / = {x e X : f{x) ^ 0}. If / is a member of £ 2 (A), then 
Supp / is defined up to a set of measure zero. 

Definition 7. Let f,g be two nonnegative functions in C 2 (X). 
Then, Supp / and Supp g are called nonoverlapping (or equiv¬ 
alently orthogonal) up to a set of measure zero if /i(Supp / n 
Supp g) = 0. 

Definition 8. A function / e C 2 {X) is said to be mixed if 

f = r-r, 

where f + = max{/, 0}, 
and f~ = max{—/, 0}, 

are called the positive and negative parts of /, respectively, and the 
following two conditions are satisfied: 

MSupp f + ) > 0 
MSupp /") > 0. 

Definition 9. An operator A is called r-potent if A r = A, 
where r is a positive integer. 

Definition 10. [9] An operator A : C 2 (Xi) —* C 2 (X 2 ) is called a 
compact linear operator (or completely continuous linear operator) 
if A is linear and if for every bounded subset M. of £ 2 (Ai), the 
image A (Ad) is relatively compact, that is, the closure A (Ad) is a 
compact subset of C 2 (X 2 ). 

Proposition 11. pi] Any two nonnegative functions in C 2 (X) 
are mutually orthogonal if and only if they have nonoverlapping 
support sets up to a set of measure zero. In other words, if /, g ^ 0 
in £ 2 (X), then (/, g) = 0 if and only if /i(Supp / n Supp g) = 0. 
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Proposition 12. If /i, / 2 ,..., f n are orthogonal functions in 
£ 2 (A), that is, /j) = 0, for all i # j, then fi, ■ ■ ■, f n are linearly 
independent. 

We next state a set of four lemmas (Lemma [T3] and [TT] are known 
in the literature (TO] while Lemma fl5l and fTfil are our contributions) 
along with the corresponding proofs that would allow us to set the 
context of the problem solved in this paper. 

Lemma 13. Let B : £ 2 (Ai) —*■ £ 2 (T 2 ) be a nonnegative oper¬ 
ator such that B/ 0 = 0 for some f Q > 0 in C 2 (Xi). Then B = 0. 

Proof. We note that 


Bf 0 = 0 

=^<B f Q , g) = 0, for all 0 in £ 2 {X 2 ) 

=></„, B*^> = 0 for all g > 0 in C 2 (X 2 ) 



=>f 0 (x)(B*g)(x ) = 0 a.e. on X\. 


However, f Q [x) > 0 a.e., and therefore 


(B *g)(x) = 0 a.e. for all g ^ 0 
B*^ = 0 for all g ^ 0 
=> B* = 0 


=> B = 0. 


□ 


Lemma 14. If the kernel of a nonnegative operator A (denoted 
by Ker(A)) contains a nonzero nonnegative function, then A is 
decomposable. 

Proof. Let h be a nonzero nonnegative function in Ker(A). Then, 
h > 0 on Supp h, implies that Ah = 0 for h > 0 in £ 2 (Supp h). 
This however, due to Lemma fl3l implies that A = 0 on £ 2 (Supp h). 
Therefore, A is decomposable. □ 

The above lemma implies that any operator A which has a 
nonzero nonnegative function in its kernel is decomposable. For 
the rest of this paper, we therefore assume that Ker(A) has no 
nonzero nonnegative functions. 
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Lemma 15. If a function / belongs to the range space R( A) 
of an r-potent operator A, then both Re(/) and Im(/) also belong 
to R( A). 

Proof. Since / e R( A) and A r_1 = I on R( A), we have 

(1) A r_1 / = / 

(2) =► A r_1 (Re(/) + z Inr(/)) = Re(/) + i Im(/) 

(3) => A r_1 (Re(/)) = Re(/) and 

(4) A r_1 (Im(/)) = Im(/), 

where Eqns. (j3j) and (Ill) follow by comparing the real and imaginary 
parts of Eqn. (J2]) as A r ~ 1 is a linear operator. However, Eqns. (131) 
and (|H) also imply that both Re(/) and Im(/) belong to R( A). □ 

Thanks to the above lemma, we shall restrict our focus in the 
rest of this paper to real functions only. 

Lemma 16. A compact r-potent operator on a Hilbert space 
has a finite-dimensional range space. 

Proof. We prove this lemma using the method of contradiction. 
Let A be a compact r-potent operator on a Hilbert space £ 2 (A) 
and let its range space R( A) be infinite-dimensional. It is easy to 
see that R( A) is a closed, and hence, complete subspace of C 2 (X). 

Now, let {e n } nE pj be an orthonormal sequence in R{ A). Then, 
A r-1 e n = e n , Vn. Also, {e n } nS N is a bounded sequence in R( A). In 
what follows, we will show that {Ae n } ne ^ does not have a conver¬ 
gent subsequence, contradicting the compactness of A. 

We start by noting that since ||e n — e m \\ = y/2 , Vn i= m, the 
sequence {e n } n6 pj, which in turn is equal to {A r ' 1 e n } nS N due to 
r-potence, cannot have any convergent subsequence. Now, sup¬ 
pose the sequence {Ae n } nE N has a convergent subsequence, say 
{Ae nA .}fc e N. Then, since R( A) is closed, {Ae nfc } ne pj must converge 
in R( A). We can therefore write 

(5) Ae nfc Aa: (say) 

(6) => A r ~ l e nk —► A r_1 a: (since A r_1 is continuous) 

(7) => e nk —»• A r_1 a: (since A r_1 = I on R( A)), 

which is a contradiction. It therefore follows that the sequence 
{Ae n } ne H does not have any convergent subsequence, which is a 
contradiction to our assertion that A is compact. □ 
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Since our focus in this paper is on nonnegative compact r-potent 
operators, we will assume throughout the paper that r-potent op¬ 
erators have finite-dimensional ranges. 

To summarise, Lemmas fl4l fl5l and fl6l imply that for the rest of 
this paper, we need to focus only on operators with (1) no nonzero 
nonnegative functions in their kernels, (2) real basis functions and 
(3) finite dimensional range spaces (say of dimension N). 

3. Characterization of Range Spaces of t-Potent 

Operators 

We show in this section that it is always possible to choose the ba¬ 
sis functions of the range space of a nonnegative compact r-potent 
operator to be all nonnegative and mutually orthogonal. The key 
word here is “nonnegative” because if we simply require an orthog¬ 
onal basis, we can obtain it using (say) Gram-Schmidt orthogonal- 
ization procedure. Typically, any orthogonal basis obtained using 
Gram-Schmidt procedure may however contain mixed functions. 
Obtaining a basis with all nonnegative yet orthogonal functions is 
completely nontrivial and is the main subject of this section. The 
rest of this section is organized as follows. 

We first present a result (stated as Lemma |T7] and Corollary fl8ll 
to show that we can always construct a basis of 72(A) comprising of 
only nonnegative functions in £ 2 (A) from any given basis of R( A). 
We then provide a set of three Lemmas (Lemma [191 [20] and I2T]) to 
demonstrate that if any pair of nonnegative basis functions have 
overlapping supports such that the measure of overlap is positive, 
then we can always obtain an alternate basis where the two ele¬ 
ments with overlapping supports are replaced by certain alternate 
basis functions with nonoverlapping supports. We finally present a 
stepwise algorithm (Theorem [22]) that uses Lemma EH to show that 
we can systematically replace all basis functions with overlapping 
supports in R( A) by alternate nonnegative basis functions whose 
supports are nonoverlapping. 

We start with Lemma [171 (generalization of [111 Lemma 3.3] from 
idempotent to r-potent case) and Corollary [T8l 

Lemma 17. If the kernel Ker(A) of a nonnegative r-potent 
operator A does not contain any nonzero nonnegative function, 
then for every mixed function / e R( A), the positive and negative 
components (namely f + and /", respectively) of / must also be 
elements of 72(A), that is, / + , /“ e 72(A). 
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Proof. Since / e R{ A), we have 

(8) A r - l f = f 

(9) - 1 (/+-/-) = /+-/- 

(10) => A r_1 / + - /+ = A"'” 1 /” - r = h (say) 

(11) => A r_1 / + = f + + h and 

(12) A r_1 /” = r + A 

Therefore, / + + h and f~ + h are both elements of R(A). However, 

(13) Ah = A(A r_1 / + - / + ) = A7 + ^ A/ + = A/ + - A f + = 0, 

that is, h is an element in Ker (A). Since there is no nonzero non¬ 
negative element in Ker (A) (recall Lemma fTTl . h must be either 
zero or a mixed function. We next show that h cannot be a mixed 
function. Specifically, we note that 

(14) / = A r -'f 

(is) =>/ + -r = A r - 1 / + -A r - 1 /" 

Since A r-1 /“ ^ 0 and f~ = 0 on Supp / + , Eqn. (1T5]) implies 

(16) A r ~ 1 f + ^ f + on Supp f + . 

Moreover, since f + = 0 on (Supp f + ) c while A r-1 / + ^ 0, we have 

(17) A r ~ 1 f + > /+ on (Supp f + ) c . 

Combining Eqn. (TT6l) and (fTTj) and noticing that 

(18) Supp f + u (Supp f + ) c = A, 
we get 

( 19 ) A^r ^ f + 

(20) => h = A r "7 + -f + > 0. 

Therefore, h cannot be a mixed function, and hence, h must be 
zero. Consequently, both f + = A r ~ 1 f + and f~ = A r ~ 1 f~ (from 
Eqns (HID and (fT21) . respectively) are elements in R( A). □ 

Corollary 18. Let {ei, e 2 ,..., e^} be a set of basis functions in 
the N-dimensional range space 77(A) of a nonnegative compact r- 
potent operator A. Then, there must exist an alternate set of basis 
functions (e), e ' 2 ,..., e' N } such that ei, Mj, are all nonnegative. 
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Proof. The proof follows in a straightforward manner by noticing 
that all the mixed functions (with both positive and negative com¬ 
ponents) among {e l5 e 2 ,. .., e N } can be replaced (due to Lemma fTTj) 
by their respective positive and negative components. By dropping 
any linearly dependent functions in the set so generated, we obtain 
the required basis {e' l5 e' 2 ,..., e^}. □ 

Corollary fl8l establishes that given any basis of an r-potent op¬ 
erator, we can always find an alternate basis where all basis func¬ 
tions are nonnegative (although they may still have overlapping 
supports). We will show in Theorem [22] (stated later in this Sec¬ 
tion) that given a basis of an r-potent operator with all nonnegative 
elements, it is further possible to construct an alternate basis where 
all basis functions are both nonnegative and mutually orthogonal 
(that is, with nonoverlapping supports). For proof of Theorem [221 
we require an important Lemma (stated later as Lemma I2TD . which 
shows that if any pair of basis functions has overlapping supports, 
then it can be replaced in the basis by alternate basis functions 
that have nonoverlapping supports. However, before proceeding to 
the proof of Lemma I2T1 we shall need the following two Lemmas as 
well. 

Lemma 19. [TOj Given two linearly independent functions 
/, g e £ 2 (fb), we can always construct a mixed function u = u + — u~ 
such that Supp u + and Supp u~ are nonoverlapping up to a set of 
measure zero. Morevover, u + ,u~ so generated would be linearly 
independent. 

Proof. Consider the following two subsets of real numbers 
/Ci = {k : f - kg > 0} 

1C 2 = {k : f ^ kg < 0}, 

and note that /Ci ^ (j) because 0 e /Ci and /C 2 4> because K 2 = 

(p => / > kg,\/k e M, which, of course, is not possible. Now, let 
ki = Supremum K,\ and k 2 = Tnfimum /C 2 . If Aq and k 2 are finite, 
then Aq ^ k 2 as then / and g would be linearly dependent, which is 
not true. Therefore, we must have Aq < k 2 (Note that since /Ci 4>, 
Aq cannot be infinity when k 2 is infinite). If we chose any real scalar 
p such that Aq < p < k 2 , then p $ /Ci because k\ < p and p $ /C 2 
because p < k 2 . Therefore, / — pg Jf- 0 and f — pg -|c 0, and hence, 
/ — pg is a mixed function. Let / — pg = u + — u~ , where u + is the 
positive component and u~ is the negative component of this mixed 
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function. Then, /i(Supp u + n Snpp u~) = 0, and consequently, 
{u + ,u~} are orthogonal and hence linearly independent. □ 

Lemma 20. Consider a mixed function u = e i — pe 2 = u + — re¬ 
generated from linearly independent functions e\, e 2 e R( A). Then, 
each of the sets {ei, u + , u~} and {e 2 ,u + ,u~} is linearly independent. 

Proof. We can divide the proof into the following cases: 

Case 1: The following three conditions hold: 

//(Supp e\ n Supp u + ) = 0, 

//(Supp e\ n Supp u~) = 0, 

//(Supp u + n Supp u~ ) = 0. 

In this case, {e\, u + }, {ei, u~}, {u + ,u~} are linearly indepen¬ 
dent. Therefore, {ei,w + ,'U - } is also linearly independent. 
Case 2: The following three conditions hold: 

yr(Supp e\ n Supp u + ) = 0, 
yu(Supp e\ n Supp u~ ) > 0, 
yu(Supp u + n Supp u~ ) = 0. 

In this case, {ei,-u + } is linearly independent. Suppose e\ = 
au~ for some a, then a > 0 and au~ — pe 2 = u + — u~. 
Therefore, 


(a + l)u — u + = pe 2 , 

OL + 1 _ 1 1 

=^e 2 = - u - u , 

p p 

which is a contradiction as e 2 ^ 0. Therefore, {ei, u ~} is lin¬ 
early independent, and consequently, {ei, u + , u~} is linearly 
independent. 

Case 3: The following three conditions hold: 

yu(Supp e\ n Supp u + ) > 0, 
yu(Supp e\ n Supp u~) = 0, 

/i(Supp u + n Supp u~) = 0. 


This case is similar to Case 2. 
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Case 4: The following three conditions hold: 

/i(Supp e\ n Supp u + ) > 0, 

/i(Snpp e\ n Supp u~ ) > 0, 

/r(Supp u + n Supp u~) = 0. 

Since //(Supp e\ n Supp u + ) > 0, we have e\ ¥= 0 but 
u~ = 0 on Supp ei nSupp u + . In addition, since /i(Supp e\ n 
Supp u ~) > 0, we have ei ^ 0 but u + = 0 on Supp u~ n 
Supp e±. Therefore, {ei ,u~} and {ei ,u + } are linearly inde¬ 
pendent. Consequently, {ei, u + ,u~} is linearly independent. 

Replacing e\ by e 2 in the above statements, we get linear indepen¬ 
dence of {e 2 ,u + ,u~} as well. □ 

We finally proceed to Lemma EH which eventually forms the 
basis of our proof for Theorem EH 

Lemma 21. Let £ = {e±, e 2 , ..., e^} be a set of nonnegative 
basis functions of R( A) such that 

(a) There exist two distinct functions e r ,e s e £ with identical 
supports, that is, Supp e r = Supp e s . Then, it is possible 
to get an alternate basis for R(A) by replacing e r ,e s in £ 
by two new functions e' r , e' s e R( A) such that Supp e! r and 
Supp e' s have at most partial overlap. 

(b) There exist two distinct functions e r , e s £ £ with Supp e r c= 
Supp e s . Then, we can construct at least two new basis 
functions e' r ,e’ s e £ with orthogonal supports. 

(c) There exist two distinct functions e r ,e s e £ with their re¬ 
spective supports satisfying the following three conditions: 
/i(Supp e r n Supp e s ) > 0, /i(Supp e r n Supp e^) > 0 and 
/i(Supp e c r n Supp e s ) > 0. Further, let e r and e s be linearly 
dependent on Supp e r n Supp e s . Then, we can construct 
three orthogonal nonnegative basis functions in R( A) from 
e r and e s . 

(d) There exist two distinct functions e r ,e s 6 £ with their re¬ 
spective supports satisfying the following three conditions: 
^(Supp e r n Supp e s ) > 0, ^i(Supp e r n Supp e^) > 0 and 
^(Supp e c r n Supp e s ) > 0. Further, let e r and e s be linearly 
independent on Supp e r n Supp e s . Then, we can construct 
four orthogonal nonnegative basis functions in R{ A) from 
e r and e s . 
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Proof. We prove the four parts of the Lemma separately: 

(a) Since the set {e r ,e s } is linearly independent, we can con¬ 

struct a mixed function u from e r and e s using Lemma [T9l 
Specifically, let e r — pe s = u + — u~, where, as previously, 
u + and u~ are the positive and negative components of u. 
Then, by Lemma ITTf we get {u + ,u~} R( A). Moreover, 

if both u + and u~ do not lie in the linear span of the re¬ 
maining basis functions £ — {e r ,e s }, then we can replace 
{e r , e s } in £ by {u + , u~} and the resulting set would be the 
required new basis. However, if either of u + or u~ lies in 
the linear span of the remaining basis functions £ — {e r ,e s }, 
then we can replace {e r , e s } in £ by the linearly independent 
set {e r , u + ,u~} (due to Lemma EUD. so that the resulting set 
£' = £ — {e s } + {u + ,u~} spans the entire R( A). However, 
£' contains N + 1 functions and therefore cannot be a ba¬ 
sis of R( A), which is iV-dimensional. Removing any of the 
functions in £' that is linearly dependent on the remaining 
functions in £' would yield the required new basis. Finally, 
note that u + and u~ cannot both be linearly dependent on 
the remaining functions £ — {e r , e s } because in this case, e s 
would be linearly dependent on £ — {e s }, contradicting the 
premise that £ is a set of basis functions. 

(b) We shall further split this part into two cases: 

Case 1: {e r ,e s } is linearly dependent on Supp e r 

Let e s = f + ae r , where / ^ 0 and Supp / n Supp e r = 
c i i. Since / = e s — ae r and both e s ,e r e R( A), we 
must have / e R( A). Therefore, {/, e r } e R( A) are 
nonnegative with orthogonal support sets and are hence 
linearly independent. Consequently, we can enter / and 
remove e s from £ to obtain a new set of basis functions. 
Please note that f £ — {e s } as otherwise e s would lie 
in the linear span of £ — {e s }, contradicting our premise 
that £ is a set of basis functions. 

Case 2: {e r ,e s } is linearly independent on Supp e r 

Let e s = f + g, where / ^ 0, g ^ 0, Supp g = Supp e r 
and Supp / = Supp e s n (Supp e r ) c . Then, 

— pe r = f + g — pe r = f + (g - pe r ) 



>o 


<o 


= u 


+ 


— u 
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where p is chosen such that g — pe r < 0 (Note that we 
can always choose such a p because otherwise g — ke r ^ 
0, for all k, which is not possible). Thus, u + = f and 
u~ = — (g — pe r ) are the positive and negative com¬ 
ponents of the mixed function e s — pe r . Consequently, 
from Lemma UTl il + = f e R( A) and u~ = — {g—pe r ) e 
R( A). However, / e R( A) implies g e i?(A) (since 
e s = f + g e R( A)). Moreover, since {e r ,g} is linearly 
independent, e r — p'g = w + — w~ is a mixed vector in 
R( A) (for some j/), and therefore, {w + , w~, /} Q R( A) 
are linearly independent functions which are nonnega¬ 
tive and orthogonal. 

(c) Let 


e r = f r + h r , where 

fr = Cr 


ISupp e r —(Supp e^nSupp e s ) 


and 


J'l/'i• c. 


Further, let 


e s fs T h s 

= f s + ah r where 
h s = ah r L „ and 

Ibupp e r nbupp e s 


fs 


Supp e 3 — (Supp e r nSupp e s )' 


Then, consider 


Cr G s f r T h v h r 

a a 


= fr 


H-ve part 


fs 


a 


-ve part 


which, due to Lemma UTl implies f r ,f s e R( A), and conse¬ 
quently, h r e R( A). Therefore, / r ,/ s ,/i r e R(A) are non¬ 
negative functions with orthogonal support sets. 
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(d) Let 


e r = fr + h r , and 
e s = fs + h a where 

fr = 6 


// j Gr 

fs 

he Gc 


Supp e r — (Supp e r nSupp e s ) 


Supp ernSupp e a 


ISupp e 3 — (Supp ernSupp e s ) 


and 


I Supp e r nSupp e s * 

Further, let a be a real scalar satisfying h r > ah s (note that 
such an a would always exist). Then, 

c r f r T h r cy^fs exhg 

fr T h> r ah s 

+ve part -ve part 

is a mixed function, and therefore {f r + h r — ah s , f s } Q 
R( A). However, f s e R(A) implies (since g s e -R(A)) that 
h s = g s — f s e R[ A). Similarly, by considering e s — ac r , 
we can show that f r ,h r e R(A). Finally, since {h r ,h s j is 
linearly independent, we can construct a mixed function x = 
h r — /3h s = x + — x~ (where x + and x~ are the positive and 
negative components, respectively, of x), for some f3, so that 
{x + ,x~, f r , f s } are four nonnegative orthogonal functions in 
R(A) which are all obtained from {e r , e s }. 

□ 


Note that we can apply the analysis of the cases l2l[Fa) andl21fb) 


respectively, to (a) any pair of basis elements in R(A) with identical 
supports, and (b) any pair of basis elements in R( A) where sup¬ 
port of one of the element is contained in the support of the other 
element, so that, eventually, every pairwise overlap of the supports 
of basis functions is partial in nature. The cases of pairwise partial 
overlap among supports of basis functions has been covered in parts 
We next state the main result of this section as 


and 121 fdl 


Theorem [22] and present an algorithm that deploys Lemmas l2i|l(c)| 
and [21[fd) in an iterative manner to systematically replace every 
pair of basis functions with partially overlapping supports in the 
basis by certain new set of basis functions that are nonnegative 
and have nonoverlapping supports, so that eventually all the basis 
functions are nonnegative and have nonoverlapping supports. 
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Theorem 22. Let £ = {ei, e 2 ,..., ejv} be a set of nonnegative 
basis functions in the range space 7Z(A) of a nonnegative compact 
r-potent operator A. Then, there must exist an alternate set of 
basis functions £' = {e^ ,e' 2 ,.... e' v } such that e'-, for all j, are 
nonnegative and have nonoverlapping supports. 

Proof. We prove this theorem by providing a constructive algo¬ 
rithm to obtain a basis with all elements nonnegative and having 
nonoverlapping supports from any given basis with all nonnegative 
elements: 


Algorithm 

Given the set of basis function £, we create two groups of func¬ 
tions: Eorth containing functions that are all nonnegative and or¬ 
thogonal, and £ NonOrth containing the remaining functions that may 
have partially overlapping supports. Naturally, the sum of the num¬ 
ber of functions in Eorth and £NonOrth would be N (that is, the 
dimension of R{ A)). The key idea behind our algorithm is to ran¬ 
domly select functions from ENonOrth , one at a time, and convert 
them into function(s) in R(A) that have nonoverlapping supports 
with all functions already in Eorth and can therefore be included in 
Eorth■ This conversion is applied in an iterative fashion so that even¬ 
tually all functions are nonnegative and have mutually orthogonal 
supports. Specifically, our algorithm consists of following steps: 

Step 1 : We begin with a single basis function, say e\, in Eorth ■ We 
initialize two counters 
Counter 1: Number of functions in Eorth 
Counter 2: Number of functions in ENonOrth 
with Counter 1 = 1 and Counter 2 = N — 1. 

Step 2 : We randomly pick one basis function, say e 2 , from ENonOrth- 
Then, Supp e 2 is either orthogonal to Supp e\ or partially 
overlaps with Supp e\. We shall consider this scenario in 
three mutually exclusive and exhaustive cases: 

Case I: e\ and, e 2 are orthogonal 

In this case, 62 is removed from ENonOrth and inserted 
in Eorth- 

Case II: {ei, 62 } is linearly dependent on Supp e± n Supp 62 
By Lemma l21jlfc) [ we can get three orthogonal nonneg¬ 
ative functions, say / 1 , / 2 , /b, in -R(A), from e\ and e 2 - 
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Therefore, we can replace e\ in Sorth by {fi, f 2 , f 3 }, and 
drop e 2 from SNonOrth■ In doing so however, the to¬ 
tal number of functions in Sorth and SNonOrth becomes 

3 + (TV — 2) = TV + 1. Since the range space 72(A) is 
TV-dimensional, one of the functions in SNonOrth must 
be linearly dependent on the remaining TV functions. 
We shall drop this function from SNonOrth so that Sorth 
and SNonOrth are together a set of basis functions. 

Case III: {ei, e 2 } is linearly independent on Supp e^Supp e 2 
By Lemma [TI Id) we can get four nonnegative orthog¬ 
onal functions, say < 71 , g 2l < 73 , g± in 72(A) from e\ and e 2 . 
Therefore, we can replace e\ in Sorth by {gi, g 2 , g 3 , g^}, 
and drop e 2 from SNonOrth■ In doing so however, the to¬ 
tal number of functions in Sorth and SNonOrth becomes 

4 + (TV — 2) = TV + 2. Since the range space 72(A) is TV- 
dimensional, there would be two functions in SNonOrth 
that are linearly dependent on the remaining TV func¬ 
tions. We shall drop these two functions from SNonOrth 
so that Sorth and SNonOrth are together a set of basis 
functions. 

Our analysis for Cases I, II and III show that Counter 1 
would now become either 2, 3 or 4, that is, Counter 1 
increases by at least one in Step 2 and, correspondingly, 
Counter 2 reduces by at least one. 

Step 3 : We select another function, say e 3 , from SNonOrth- Then, 
e 3 would be either orthogonal to all functions in Sorth or 
Supp e 3 would partially overlap with one or more of the 
functions in Sorth■ If e 3 is orthogonal to all functions in 
Sorth, it can be directly removed from SNonOrth and inserted 
in Sorth, thereby increasing Counter 1 by 1 and reducing 
Counter 2 by 1. If, on the other hand, Supp e 3 overlaps 
with one of the functions in Sorth, say h, then we shall have 
the following two cases: 

Case I: {h, e 3 } is linearly dependent on Supp h n Supp e 3 
By Lemma |5](c) we can get three orthogonal basis 


functions from h and e 3 , namely, /(, / 2 , / 3 such that 


Supp f[ c: Supp h n (Supp e 3 ) c 
Supp f 2 <= Supp e 3 n (Supp h) c 
Supp /g ci Supp e 3 n Supp h, 
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respectively. 

Case II: {h, e 3 } is linearly independent on Supp hr)Supp e 3 
As {h, e 3 } is linearly independent on Supp h n Supp e 3 , 
by Lemma 121 Id) we would get four orthogonal basis 


functions, namely, g [, g ' 2 , g 3 , g' A such that 


Supp g\ <= Supp h n (Supp e 3 ) c 
Supp ^ — Supp e 3 n (Supp h) c 
Supp g' 3 Q Supp e 3 n Supp/i 
Supp g' A ci Supp e 3 n Supph, 


respectively. 

If Supp e 3 overlaps with another function (say h') in Sorth, 
then it would be on Supp e 3 n (Supp/?.) c (recall {h, h'} have 
orthogonal supports). In addition, in Case I above, (Supp hi 
nSupp e 3 ) cc Supp f 2 and therefore we can repeat the same 
argument as above for {h r , f 2 } rather than {h. e 3 } to obtain 
three new orthogonal functions, and so on. In Case II above, 
We have (Supp h! n Supp e 3 ) Supp g 2 and therefore we 
can repeat the same argument as above for {h ', g 2 } rather 
than {h, e 3 } to obtain four new orthogonal functions, and 
so on. The above process is repeated till all the overlaps 
between e 3 and functions of Sorth are eliminated. The re¬ 
sulting orthogonal functions are aggregated in the new Sorth 
and any set of functions in £NonOrth that are all linearly de¬ 
pendent on the remaining N functions in Sorth and SNonOrth, 
are dropped. 

By performing the above step, we get a new Sorth where number 
of functions in definitely greater than the number of functions in 
£orth obtained after Step 2. That is, due to Step 3, Counter 1 
increases by at least one, Counter 2 decreases by the same number 
as the increase in Counter 1 and the functions in Sorth u SNonOrth 
still constitute a basis. 

Finally, note that repeating Step 3 above for every remaining 
function in SNonOrth, one at a time, we will eventually get Counter 
1 = N (recall that number of basis functions cannot exceed the di¬ 
mension of R( A), which is N ). The Sorth so obtained would be the 
required basis of R( A) with all functions nonnegative and mutually 
orthogonal. Therefore, the proof for Theorem (22} is complete. □ 
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4. DECOMPOSABILITY OF r-POTENT OPERATORS 

Our main result in this section is stated as the following theorem: 

Theorem 23. If there exists a basis {ei, e 2 ,..., e^} in the 
range space R( A) of a nonnegative compact r-potent operator A 
with r ^ N such that e 3 , for all j, are nonnegative and have non- 
overlapping supports, then A must be decomposable over some 
support set U. 

Proof. We will first prove the above theorem for the case where 
r = 3 and will then generalize the proof to r > 3 case. 


Case: r = 3 


We start by noting that Ae* e R( A), for i = 1,..., N. Therefore, 


( 21 ) 


Ae, — a^ei + 0(262 + • • • + Q(jyejV; 


for some a\,..., a l N 5= 0 (since A is a nonnegative operator and Cj, 
for all j, are nonnegative orthogonal). 

We next claim that Ae, must lie in the linear span of exactly one 
e 3 {j = 1 ,... ,N), that is, there must exist only one j 0 such that 
a l j =£ 0, for j = j 0 , and ck*- = 0, for all j ^ j Q . We can prove this 
claim as follows: 

First, consider that 

(22) 6j = A ei = oyAe^ + OL 2 Ae 2 T ■ • • + cvjyAcjy. 

Now suppose that two coefficients, say a* and o®, are both greater 
than zero, and rewrite Eqn. (1221) as 



j^=P,Q 

However, since both Ae p , Ae g e R( A), we can write 


(24) Ae p — afei + o 262 + • • • + o^e^v 

(25) Ae q = a\e 1 + a\e2 + ■ ■ • + a q N e N . 


Substituting Eqns. (041) and (1251) into Eqn. (125j) . we get 

ei = o® (ccfei + a p 2 e 2 + ■ ■ • + u p N e N ) 

+ of q (( x^e 1 + 0(262 + ■ ■ • + o ^-6 jv ) 





(26) 
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As {ei, e 2 ,..., ejv} is linearly independent and a* are all nonnega¬ 
tive, we must have 

a i = a 2 = ''' = a i-i = a i +1 = • • • = oi p N = 0 and 

( 27 ) a\ = a% = ■ ■ ■ = ol\_ x = a g i+1 = ■ ■ ■ = a q N = 0 

which implies 

Ae p = a p ei and 

( 28 ) Ae q = a g ei 

where both a p ,aj > 0 (otherwise, A would not be 3-potent). 
Therefore, 


(29) e p = A 2 e p 

(30) = of A ei 

(31) = QL P {oi\^& i + OL2&2 + ' ' ' + 



which, again using linearly independence of {ei, e 2 ,..., e^} implies 
a p a q = 0 . ffowever, this is a contradiction to our assumption that 
both a p , of > 0 , and therefore, our claim that Ae* must lie in the 
linear span of exactly one ej is proved. In other words, we must 
have a unique j for which 

(33) Aej = ci-jCj, where ct* > 0. 

Note that if dim R(A) = 2, then {ei, e 2 } is a nonnegative orthog¬ 
onal basis of /2(A). Therefore, due to the above proved claim, Aei 
is equal to either ae\ or a'e 2 , for some a, a' > 0. Similarly, Ae 2 is 
equal to either /5e 2 or j3'e i, for some /3,/T > 0. 

But 

(34) Aei = ol€\ 

(35) => ei = A 2 ei = q.A&i = a 2 ei 

(36) => a = 1 (since a > 0), 

which yields 

(37) 


Ae! = e\. 



DECOMPOSABILITY OF NONNEGATIVE r-POTENT OPERATORS ON C 2 {X§ 


In addition, Aei = ae i also implies Ae 2 = f3e 2 (since Ae 2 = j3'ei 
would imply e 2 = A 2 e 2 = (3 'Aei = a/3'e i, which would be a con¬ 
tradiction), which in turn implies 

(38) e 2 = A 2 e 2 = /3Ae 2 = f3 2 e 2 , 

which, along with the condition (3 > 0, implies Ae 2 = e 2 . That is, 
A = I on R( A). In other words, A is a 2-potent (or idempotent) 
operator which is already known to be decomposable. Therefore, 
we may assume that Aei ¥= aei and Ae 2 ^ (3e 2 . 

Consequently, 

(39) Aei = a'e 2 

(40) => ei = A 2 ei = a'Ae 2 = a'/3'e i, and 

(41) Ae 2 = (3'e 1 

(42) => e 2 = A 2 e 2 = (3 'Aei = (3'a'e 2 . 

Therefore, applying A to ei yields e 2 and applying A to e 2 , in turn, 
yields ei. This implies that if we need decomposability, we must 
have dim R( A) > 2. 

In case dim R( A) > 2, let us consider ei + Aei. Then, 

(43) //(Supp (ei + Aei)) > 0. 


u 


Also, since dim R( A) > 2 , 

(44) 

Further, 

(45) 

(46) 

(47) 

(48) 


n(U c ) > 0 . 

(ei + Aei, Xu c ) = 0 

=> (Aei + A 2 ei, xu c ) = 0 
=► <ei + Aei, A *xu°) = 0 


(ei + Aei) (a;) A*x^ c (x)fi(dx) = 0 


x 


However, since (ei + Aei)(x) > 0 on U , Eqn (1481) implies that 

(49) A *Xu c (x) = 0 a.e. on U 

(50) ^<A* Xw c, Xw > = 0 

(51) =* (Xu*, A Xw > = 0 , 

and therefore, as per Defn. [4l the 3-potent operator A is decom¬ 
posable. We next generalise the above proof to the r > 3 case. 
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Case: r > 3 

We start by noting that Ae* e R(A), for i = 1, 2,..., N, and 
therefore 

(52) Ae, = a[ei + 0^262 T ■ • • + cu^ejv, 

for some a\,... , a l N ^ 0, where at least one of a \,..., a l N is strictly 
greater than zero or else A would not be an r-potent operator. 

We next claim that Ae* lies in the linear span of exactly one 
ej (j = 1 , 2 ,..., 1 V), that is, there exists a unique j„ such that 
ctj > 0, for j = j 0 and cc* = 0, for all j ^ j„. We prove this claim 
by the method of contradiction. We assume that there exists two 
coefficients a l p and a l q such that both a l p , a l q > 0. Then, 

(53) A e-i = a l p e p + a q e q + ^ ot^j 
and therefore 

A 2 ej = a p Ae p + a l q Ae q + ^ a* A e-j 

(54) = oj p {pie' p + NNOT) + a^e'g + NNOT) + NNOT, 

where e p and e' q are some basis functions whose coefficients are 
nonzero in the expansions of Ae p and Ae q , respectively (recall that 
at least one of the coefficients must be strictly positive in every 
expansion else A would not be r-potent), pi,qi are those strictly 
positive coefficients, and NNOT (NonNegative Other Terms) is a 
general term introduced for the sake of convenience to denote any 
sum of nonnegative functions whose exact value is not relevant for 
further analysis in this proof. 

In a manner similar to expansions of Ae* and A 2 e,; in Eqns. (HuH) 
and §H), respectively, we can write expansions of A 3 ej, • • • , A r_2 ej 
as 

A i e i = a p A 2 e p + a l q A 2 e q + ^ a)A 2 ej 

= a p (piAe' p + NNOT) + a*(giAe^ + NNOT) + NNOT 
= c4(pi(;p 2 e" + NNOT) + NNOT) 

+ ctq{qi(q 2 e" + NNOT) + NNOT) + NNOT 


(55) 
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(56) A r ~ 2 ej = aj, A r ~ 3 e p + a\A r ~ 3 e q + NNOT 

(57) = OL p p x p 2 ■ . . p r - 3 e p + a l q qiq 2 ... q r - 3 e q + NNOT 

where a p ,p\,p 2 - ; . . . ,p r -3 and a q , qi, q 2 ,..., q r -3 are all chosen to 
be positive (note that we can do so because for any ej, A ej must 
have at least one coefficient a 1 * > 0 else A would not be r-potent). 
Applying A to Eqn. (1571) yields 

(58) ei = A!'~ l e-i 

(59) = a^p^ ... Pr- 3 Ae p + a\q x q 2 ... q r - 3 Ae q + NNOT 

In Eqn (1591) . at least one of the coefficients in expansion of both 
Ae p and Ae q is nonzero. Without loss of generality, let the term 
corresponding to the nonzero coefficient of Ae p be /3 m e m and the 
term corresponding to nonzero coefficient of Ae q by 7 n e n . Then, 

(60) ei = {a l 0 p\p 2 ... Prsl3 m )e m + (a*gi ? 2 • • • 9n-37n)e n + NNOT. 

As the coefficients of both e m and e n in the above equation are 
strictly greater than zero, we must have both m = i and n = i 
(otherwise, due to linear independence of e x ,... ,e^, all the coeffi¬ 
cients would be zero). Therefore, 

(61) Ae p fii&i 

since Ae p will not involve any other nonzero term because then, as 
earlier, there will be a contradiction to our assertion that ej for all 
j, are linearly independent. By a similar argument as above, we 
can show that 

(62) Ae q = 7 i6i. 

Since Eqns. (1561) and (1571) sets 

(63) A r ~ 3 e p = pip 2 .. ,p r _ 3 e p + NNOT, and 

(64) A r ~ 3 e q = q x q 2 ... q r - 3 e q + NNOT, 

we can now write 

A r ~ 2 e p = pip 2 .. .p r _ 3 Ae p + NNOT 

(65) = P 1 P 2 ■ ■ -Pr-3l3iei + NNOT and 
A r ~ 2 e q = gi q 2 . . . q r - 3 Ae q + NNOT 

= q±q 2 • • • qr-zli^i + NNOT, 


(66) 
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which further implies 


(67) e p = A r_1 e p = pip 2 • •-Pr-sAAe* + NNOT and 

(68) e q = A r-1 e g = q x q 2 ... g r _ 3 y^Ae* + NNOT. 
Expanding Eqn. ()67l) for Ae* using Eqn. ((53]), we get 



which is a contradiction to the linear independence of e p and e q as 
P 1 P 2 ■ ■ -PrsPiCXq > 0. Therefore, the following must hold 

(70) Ae* = a)ej, 

for a unique j such that a* > 0. This proves our claim that Ae^ 
must lie in the linear span of exactly one e r 

Since the claim holds for all e*, it follows that Ae* lies in the linear 
span of e.j , A 2 e* = a 1 - Aej lies in the linear span of (for some k), 
and so on. However, since A r_1 ej = e*, if we keep successively 
applying A to a basis function e i) e* must repeat in a cyclic fashion 
(at least once for every r — 1 applications of A). Please note that 
for at least one e^, the frequency of this cyclic repetition should be 
r — 1, else A would not be r-potent. 

Now, if dim R( A) = r — 1, then {ei, e 2 ,..., e T _i} would be a 
nonnegative orthogonal basis for R( A) and A would not be de¬ 
composable. Therefore, in order to have decomposability of an 
r-potent (but not fc-potent for k < r) operator, the dim R( A) must 
be greater than r — 1. 

If dim R( A) is indeed greater than r — 1, then 

C 2 ^Supp (e* + Ae* + A 2 e* + • • • + A r-2 ej)^ 

would be the required decomposing space with p(U) ■ p(U c ) > 0. 
Specifically, 

(71) (ei + Ae.j + A^e* + • ■ ■ + A' 2 ej, Xu^ = 0 

(72) => (A' 1 e, + A ei + A J e, + • ■ ■ + A' 2 ej, Xu c ') = 0 

(73) => (ei + Ae* + A“e.j + • ■ ■ + A' 2 e*, A *Xu c ') = 0 

(74) => f (ei + Ae* 4-h A r ~ 2 e, : ) (x)A*xw(x)p(dx) = 0 

Jx 
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which implies 

(75) A *Xu c (x) = 0 a.e. on U 

(76) ^< A * Xwc , Xw > = 0 

(77) => <Xw=, Axw> = 0, 

which, by Defn. SI implies that the nonnegative r-potent operator 
A is decomposable. □ 
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